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Sequence (DS)- Code Division Multiple Access (CDMA). Assuming perfect estimates of channel and 
interference covariance, the receiver selects the signature that maximizes signal-to-interference plus 
noise ratio (SINR) for a desired user from a signature codebook. The codebook index corresponding 
to the optimal signature is then relayed to the user with finite number of bits via a feedback channel. 
Previously, we showed that a Random Vector Quantization (RVQ) codebook, which contains independent 
isotropically distributed vectors, is optimal (i.e., maximizes SINR) in a large system limit in which 
number of interfering users, processing gain, and feedback bits tend to infinity with fixed ratios. Dai et 
al. have analyzed the large system SINR for a matched filter with nonfading channel. Here we extend the 
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I. Introduction 

A user performance in Direct Sequence (DS)- Code Division Multiple Access (CDMA) 
depends on a signature code, which can be optimized to increase signal-to-interference plus 
noise ratio (SENR). Several works in the literature [l]-[7] have investigated a joint transmitter- 
receiver signature optimization and showed that a performance difference between optimized and 
random signatures can be substantial. However, adapting signature adds more complexity and 
requires knowledge of channel and interference covariance at both the transmitter and receiver. 
All aforementioned works assume that perfect estimates of channel and interference covariance 
are available. This assumption, especially at the transmitter, is not practical. 

A receiver typically estimates channel coefficients and interference covariance from pilot 
signals during training period. The accuracy of the estimation increases with amount of available 
pilots. The transmitter, on the other hand, is usually unable to directly estimate the forward 
channel. However, channel information may be obtained from the receiver via a feedback channel. 
Thus, accuracy of channel information at the transmitter depends on available feedback rate, 
which is normally low. In recent years, many researchers [8]-[15] have proposed feedback 
schemes in which the receiver computes and quantizes the optimal signature and relays the 
quantized coefficients to the transmitter via a low-rate feedback channel. References [10]-[15] 
consider multiantenna systems where spatial signatures are optimized and quantized. Here our 
interest is signature quantization in DS-CDMA and its performance, which depends largely on 
quantization codebook and available feedback rate. 

The signature codebook is known a priori at both the transmitter and receiver. With B 
feedback bits, the receiver selects the signature vector, which maximizes the instantaneous SINK, 
from 2 ^-signature codebook and relays the corresponding index to the transmitter via an error- 
free feedback channel. References [8], [16] proposed a Random Vector Quantization (RVQ) 
codebook, which consists of independent isotropically distributed vectors and showed that the 
RVQ codebook is optimal (i.e., maximize the SINR over all codebooks) in a large system limit 
in which number of users K, processing gain N, and feedback bits B tend to infinity with 
fixed K = K/N and B = B/N. The upper bound on asymptotic SINR for single-user matched 
filter was derived in [8]. Reference [8] also considered a minimum mean square error (MMSE) 
receiver and derived an approximation for a large system SINR. The large system performance 
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was shown to predict the performance of a finite-size system well for small B. 

Recently, [9] derives the exact expression of a large system SINR for RVQ with a matched 
filter and ideal nonfading channel. (Similar results for the performance of RVQ in multiantenna 
system were derived in [12].) Here we extend the results for matched filter to multipath fading 
channel and arbitrary transmit power across users. We apply similar techniques used in [9], [12] 
to derive expressions for asymptotic SINR with linear MMSE receiver. For an MMSE receiver, 
we first consider a nonfading channel and derive an exact expression for a large system SINR, 
which is a function of K and B. We remark that the expressions for MMSE receiver are not 
trivial extensions of [9], [12]. Comparison between the large system SINR and the approximation 
derived in [8], which over-estimates the performance for large B, is shown. Numerical examples 
show that the large system results estimate the performance of the finite-size system well. From 
examples shown, one feedback bit per signature element achieves close to the performance with 
unlimited feedback. 

II. System Model 

We consider a discrete-time reverse-link synchronous DS-CDMA in which there are K users 
and processing gain N. The N x 1 received vector is given by 

K 

r = ^ y/A^HkSkbk + n (1) 

k=l 

where ^/A^ is the amplitude of user k, Hk is the N x N channel matrix for user k, Sk is 
the X 1 signature vector for user k, hk is the transmitted symbol for user k, and n is the 
additive white Gaussian noise with zero mean and covariance cr^/. For ideal nonfading channel, 
Hk — I. For frequency-selective channel, we assume that the symbol duration is much longer 
than the delay spread and, thus, we discard any inter-symbol interference. Assuming that each 



August 17, 2009 



DRAFT 



4 



user traverses L Rayleigh fading paths, we have 
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where fading gains for user k, hk,i, ■ ■ ■ , hk,L, are independent complex Gaussian random variables 
with zero mean and variances E\hk,i\'^, . . . , E\hk,i\'^, respectively. For flat fading channel (L — 1), 
Hk — hk,il. 

The receiver applies linear filter on the received signal to obtain the received symbol. We 
consider both matched filter and linear MMSE receivers and assume, without loss of generality, 
that user 1 is the user of interest. The matched filter for user 1 is given by 



Ci = Si 



(3) 



where we let Sk = H^Sk, which is the effective signature for user k. The matched filter is 
simple and can be a performance benchmark for a more complex receiver. The associated SINR 
is given by 



7 = 



c\RiCi 

where the interference-plus-noise covariance 



s\hIRiHiSi 



R, = E[rirl] 



where expectation is over transmitted symbols and noise and 

K 

n^r- s/AHiSibi = ^ s/Ai^HkSkbk + n. 

k=2 



(4) 



(5) 



(6) 



Assuming that bk's are independent and identically distributed (Ltd.) with zero mean and unit 
variance, we have 

K 

Ri^Yl ^^^^4 + = .^i^i-^I + all, (7) 



k=2 
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where Si is the N x {K — 1) effective signature matrix whose columns consist of s^, V/c 7^ 1 
and Ai is the (K — l)x (K — 1) where diagonal matrix whose diagonal entries are A2, . . . , Ak. 
On the other hand, a linear MMSE receive filter for user 1 is given by 

ci = R-^Si (8) 

where the received covariance 

K 

R = E[rr''] = J2 Akhsl + all, (9) 

k=l 

assuming that bk's are (i.i.d.) with zero mean and unit variance. Similar to the matched filter, 
we can compute the SINR for user 1 given by 

(3 = Ais\hIr^^HiSi (10) 

where matrix inversion lemma was used to simplify the expression. A linear MMSE receiver is 
shown to be robust in suppressing multiple- access interference [17]. We note that, for given Ri 
and Hi, the SINR for user 1 is a function of the signature si for both receivers. 

The receiver, which is assumed to have a perfect estimate of the interference covariance -Ri 
and channel matrix Hi, can optimize the signature for the desired user to maximize the received 
SINR. Ideally, the receiver sends the optimal signature back to user 1 via a feedback channel 
and the user changes the signature, accordingly. Practically, a feedback channel has limited rate 
and thus, the receiver can only relay finite number of feedback bits to the user. (We assume 
that the feedback does not incur any errors.) With B bits, the receiver selects the signature 
from a signature set or codebook containing 2^ signatures. This codebook is designed a priori, 
and is known at both the user and receiver. The performance of the optimized user depends 
on the codebook. Several work [8], [10]-[14], [18] focused on codebook design and analyzed 
the associated performance. (All except [8] are in context of spatial signature in a multiantenna 
channel.) In this work, we analyze the performance of a Random Vector Quantization (RVQ) 
codebook proposed by [8], [16]. RVQ codebook 

V^{Vi,...,V,b} (11) 

in which v/s are independent isotropically distributed with unit norm (\\vj\\ = 1). In other words, 
signature vectors in RVQ codebook are uniformly distributed on a surface of an A^-dimensional 
unit sphere. In [8], [9], [12], RVQ was shown to maximize SINR over all quantization codebooks 
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in a large system limit to be defined. Although RVQ is optimal in a large system limit [8], it 
was shown to perform close to the optimal codebook designed for a finite-size system [19]. 
Given the codebook V, the receiver selects 

Si = argmax SINR(t),). (12) 

The index of the optimal signature vector is relayed to user 1 via a feedback channel. We are 
interested in analyzing the corresponding SINR, which is a function of available feedback bits, 
for both matched filter and MMSE receiver. 

III. Large System Performance 

A. Matched filter 

We first consider ideal nonfading channel (Hi = I for all i) for which the optimal signature 
that maximizes SINR also minimizes the interference. Given the RVQ codebook V, the optimal 
signature is given by 

Si = argmin ,) = vIRivA (13) 

where / is the instantaneous interference power. Since v/s in RVQ codebook are i.i.d., the 
corresponding /(i'j)'s are also i.i.d. and thus, the associated interference averaged over codebook 
is given by 

PCX) 

Ev[/K),...,/(«2s)|i?i] = 2^ / x[GnR,ix)f-'gjin,ix)dx (14) 

Jo 

where Gj\ji^{-) and gi\R^{-) are cumulative distribution function (cdf) and probability density 
function (pdf) for I{vj), respectively. It is difficult to evaluate (O for any finite N, K, and B. 
It was shown that the interference power converges to a deterministic value in a large system 
limit in which K, N, and B all tend to infinity with fixed normalized load K = K/N and 
normalized feedback bits B = B/N [8], [9]. Applying theory of extreme order statistics [20] 
similar to [8], the large system interference power with fading channel is given by 

I^. = lim Gj' (2-^) (15) 

where the empirical eigenvalue distribution of -Ri converges almost surely to a nonrandom limit 
as (A^, K) oo with fixed K/N. Rearranging ([T5]) gives 

lim [GriRAzp=2-^. (16) 

{N,K,B)-*oo 
^ '-'rvq 
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References [9], [12] show that 



lim [Gj\i,,iz)]Tr = exp{-vl/(pM-)} (17) 



i 1° 
^rvq 



where 



^(P,Cq) = / log(l + p(A-Cq))(7K.(A)dA, (18) 

p* = arg max ^(p,Cq), (19) 



rvq ^min 



and gR^{-) is an asymptotic eigenvalue density for Ri. Equating (fT6l) and (flTI) . we have that 

*(p*,Cq)=51og(2). (20) 



/r^q satisfies 



With ©, Eq. ([H]) becomes 

*(P,/rvq) = / log(l + p(A + a2 -/-))^^^^^^t(A)dA. (21) 
= log(l + pal - pl^^) + J log(l + ^X)9s,^^sl ^^^^ 
= log(l + pal - pl^^) + z^^^^^^t (0 (23) 

where 



and «^s^A^st(-) is the Shannon transform for asymptotic eigenvalue distribution for 
Reference [21] defines the Shannon transform for a density function /x(-) as follows 

i^xh) = J log(l + jx)fxix) dx. (25) 

Suppose Sk, 2 < k < K, has independent complex Gaussian entries with zero mean and variance 
{\\sk\\ 1). The eigenvalue distribution for SiSl converges to a deterministic function 
as N,K ^ oo with fixed K [22] and we assume that empirical distribution of A2, . . . , Ak 
converges to a limit. It is shown that [21] 
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where VsiAiS^^'^ ry-transforai for the asymptotic eigenvalue distribution for SiAiSl and 

the //-transform for a distribution for random variable X is defined in [21] as follows 

Vx{l)= [ -^fx{x)dx. (27) 

With earlier assumption on the distribution for Si, [21] shows that VsiAisI^-'^^ solution to 

the following fixed point equation 

K = P^^^-' (28) 



Combining ^ UB, and (EH]), we have the first result. 

Theorem 1: The large system interference power at the output of single-user matched 
filter satisfies the following equation 

max {log(l + pal - pl^^) + KuAM^iO) " log 6(0 + 6(0 - 1} = 51og(2) (29) 

0<P<7oc>A 

^rvq ^mm 

where 6(x) is the solution to the following fixed point equation 

K = ^-Q(^) (30) 
1 - ?7Ai(a;6(x)) 

and ^ is given by (|24l) . 

Solving for I^^ requires numerical solution in most cases. However, for equal power allocation 
(Ai = A2 = ■ ■ ■ = Ak), the explicit expression for I^,^ was obtained by [9] as follows. 
Corollary 1 ( [9]): Let B* = ^"^^Igff ^'"^ for K > 1. For K > 1 and B > B*, 

^rvq = + (1 - V^)' + V^(l - ^)'-''exp(-v^- 51og(2)). (31) 
Otherwise, 1^^^ = Q + where Q satisfies the following equation 

Q = K^iQ-K)/K^-B/K_ (32) 

Thus, interference power decreases exponentially with the normalized feedback bits and is near 
a single-user performance with only few feedback bits per processing gain. The associated SINR 
for user 1 in a large system limit is then given by 

-"rvq 
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If the interfering signatures are orthogonal = /) with equal power allocation, the 

eigenvalue distribution is given by for K < 1 

9s,A,slW = ^'^(^ - 1) + (1 - K)S{X). (34) 

Evaluating (|20|) with the distribution in (|34|) . we obtain the following result. 

Theorem 2: The large system interference power I^^ for orthogonal interfering signatures 
with equal power allocation satisfies the following fixed-point equation 

(I^, - + - Cq)'-^ = ^^(1 " i^)^~^2-^ (35) 

for < < 1. 

For K ^ 1, we obtain the following approximation 

I^^^al + K2-^. (36) 

Here we see clearly that the interference power decreases exponentially with the normalized 
feedback bits. 

For a flat-fading channel (L = 1), we can combine the channel gain |/ifc,ip for user k 
with its transmit power A^. That is the diagonal matrix Ai = diag{|/i2,iP^2, • • • , \hK,i\'^^K} 
whose empirical distribution converges to a nonrandom limit. With an asymptotic distribution 
for diagonal entries of Ai, we can apply (|20l) - (|30l) to obtain the output SINR 7^^q. 

For frequency-selective fading, the signal of each users is assumed to propagate L discrete 
chip-spaced paths. The channel matrix for user k is shown in (|2l). First, we assume that L is 
finite and does not grow with A^. Thus, the number of paths per processing gain L = L/N 0. 
To compute I^^, we require the asymptotic eigenvalue distribution of -Ri ([7]). Reference [23] 
showed that the asymptotic eigenvalue distribution of Ri with L-path channels Q equals 
that of -Ri with flat-fading channels and Ai = diag{y42(X]^i \h2,i\^), ■ ■ ■ y^KC^i=i l^/r./P)}- 
Thus, a multipath interferer is asymptotically equivalent to a single-path interferer with com- 
bined gain of J2f=i \hk,i\^- For ^ oo with fixed L/N, the same result applies as long as 
lim^^^oo Ef=i E\hk,i\^ < oo, for aU k. 

The associated SINR at the output of matched filter with fading channel is given by 

7- = ^ (37) 

'rvq Too ^^'^ 
-"rvq 

where J2i=i ~^ ^^^d J^^q is obtained by Theorem [T] 
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B. Linear MMSE Receiver 

The SINR with the optimal signature averaged over the RVQ codebook is given by 

Ev[mBJL{(3{vi),...,l3{v2B)}\RuHi] = 2'' / x[Fp\R,^HAx)f'^fm^,HAx)^^ (38) 

Jo 

where fi3\R^,Hii-) and -F)3|i?i,fji(-) be pdf and cdf for the output SINR l3{vj), respectively. Similar 
to the matched filter, computing (|38l) for finite parameters are difficult. Taking the large system 
limit as A^, K, B ^ oo with fixed ratios, the SINR converges to a deterministic value 

13^^= lim Ev[max{/3i,...,/32s}|i?i,ifi] (39) 
= lim F-i, „ (1 - 2-^), (40) 

which can be shown by applying theory of extreme order statistics [20]. Reference [8] derived 
the approximation for by approximating cdf for P{vj) to be Gaussian. The approximation 
is a function of K, B, and cr^ and is good for small B. For large B, it over-estimates the actual 
performance. In this section, we derive exact expressions for /S^^q. 

We first consider the ideal channel (Hk = I, V/c). We rearrange (l40l) to obtain 

^^lim [1-F,|H,(^#=2-^. (41) 

(N,K,B)^oo 
^ f-'rvq 

Similar to [9], [12], it can be shown that 

lim [1 - Fp^n,{zp = exp{-$(p*,/3~ )} (42) 



where 



<^'(P,/?rvq) = / log(l + p(/3rvq--^))/s,A,stM dr, (43) 
•J ' n 

p* = arg max ^{p^f^^n)^ (44) 



f-'max '-'r\ 

t 



•^SiAist(') is the asymptotic eigenvalue density for SiAiS\, Si is the x (K — 1) signature 
matrix whose columns are S2, . . . , s^, and /9^ax i^ the asymptotic maximum eigenvalue of AiR^^ 
and corresponds to the SINR with infinite feedback (B oo). 

Combining (|4T]) and (|42|) . /^^'^q satisfies the following fixed-point equation 

$(p*,/3,^q)=51og(2). (45) 
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To evaluate ^{p*,/3^q), we rewrite (|43] ) as follows 

<^>(P,/?rvq) (46) 
= [ log(l + Cr)4^^^st(r) dr- [ log(l + -ir)4^^^^t(r) dr + log(l + p(/5,^q - ^)) 

(47) 

= ^5. A,5l (0 - ^S^MSI ( V^^) + log(l + p(/3- - ^)) (48) 

where 



300 



^ A ^ + ^^^vq (49) 



and z^5^^^5t(-) is the Shannon transform for asymptodc eigenvalue distribution for 
With similar steps used to derive Theorem [H we obtain the following theorem. 
Theorem 3: For B, the large system SINR is given by 



max {log(l +p(/5,^q - ^)) + Kiy^ACm) " Ki^aA^-'QK')) 



0<p< i^oii "Lflgo 



71 



(50) 



f-'max /-^rvq 

-log(e(C)) + log(e(a-2)) + 6(0 - e{a-')} = 51og(2) 
where ^ and 6(s) are given by (|49l) and (|30l) . respectively. 

For an equal-power (Ai = ^2 = . . . = A^) system, we can simplify expression for {3^^ as 
follows. 

Corollary 2: We assume ideal channels, i.i.d. interfering signatures, and equal transmitted 
power across users. For K < 1, satisfies the following equation 

log(T3|^ - ^) + (1 - ^) log(^) 



where 



wi^l) 1 - v{al) 

- v{p) + v{al) = B \og{2) (51) 



w{x) = ^{l + K + x+\/{l + K + xy-4K) (52) 



^;(a;) = + K + -y/(l + /s: + x)2 -4/^) (53) 
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and 



I - goo 2 

P = T>7^ . n T - 7^ ^ K- (54) 



For K > 1 and B < B*, satisfies the following equation 



- 1) log(^44) - ^(P) + = B log(2) (55) 



where 



K\og{^K-l) - ^K-\og{w{ai)) 



+ {K-l)\og(l-'^) + v{al)). (56) 



For > 1 and B > B*, 

/^rvq = /5max(l - 2-^[exp{iinog(ir) 



(K - 1) log( ^ , 5 - log(^(a^)) 



+ v{al)-^/R}]). (57) 



The derivation is shown in the Appendix. 

We can also derive the SINR when the interfering signatures are orthogonal. Substituting the 
corresponding eigenvalue distribution (|34l) into (|47] ) and simplifying give the following result. 

Theorem 4: For orthogonal set of interfering signature with < A' < 1 and equal power 
allocation, the large system SINR (3^^ is the solution of the following fixed-point equation 



{A.-p-alfW- - {A-p-al)r^ = (^^i^^V'''2-^. (58) 

For a system with heavy load (K ^ 1), we have 



» ^ - j;(TT^2-" (59) 

The first term on the right-hand side of (|59| ) is the single-user performance and thus, the 
performance with the MMSE receiver also increases exponentially with B, which is the same 
as the that of matched filter, which is a much simpler receiver. 
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Similar to the matched filter, a multipath fading is asymptotically equivalent to a single fading 
path with combined gain of Yl!i=i \ hk,i\^- Combining the fading gain with transmitted power, we 
obtain a new Ai. We assume that the distribution for a diagonal elements of Ai converges to a 
deterministic function given that < oo and X^Ili l^fe./P < ^ k. 

IV. Numerical Results 

Fig. [H shows the asymptotic SINR for MMSE receiver in Corollary [2] versus normalized 
feedback bit B with different normalized loads K = 0.25,0.5, 1, 1.25. As expected, the SINR 
increases with normalized feedback and decreases with normalized load. For K = 0.25, RVQ 
achieves close to the single-user performance with approximately B = 0.5 (0.5 bits per process- 
ing gain or degree of freedom). As number of interfering users increases, amount of feedback 
required also increases to achieve a target SINR. For example, B = 3 is needed for system 
with ^ = 1 to achieve close to the single-user performance. We also compare the asymptotic 
results with simulation results marked by pluses in Fig. [B We note that the large system results 
predict the performance of finite-size systems (A^ = 12) well. As N increases, the gap between 
the simulation and analytical results is expected to be closing. RVQ codebook requires an 
exhaustive search to locate the optimal signature. The search complexity increases exponentially 
with feedback bits B. (For B = 3, number of entries in RVQ codebook is 2^^^.) Thus, we do not 
have simulation results for a large B. 

Fig. [2] shows the large system performance of both receivers, which is obtained from Theo- 
rems [T]|4l with different distributions for interfering signatures. Interfering signatures either have 
i.i.d. Gaussian elements or are orthogonal. From the figure, the system with independent Gaussian 
signatures performs a bit better than that with orthogonal signatures for both linear receivers. The 
difference is more pronounced in the matched filter. Which distribution for interfering signatures 
gives the maximum performance is an interesting open problem. 

In Fig. [3l we compare the asymptotic SINR for MMSE receiver in Corollary [2] with the 
approximation derived in [8] for K = 0.75 and SNR = 10 dB. Also shown is the simulation 
results with N = 12. The large system SINR is closer to the simulated performance than the 
approximation. We also show the RVQ performance of a matched filter in Corollary \T\ [9] with 
that of MMSE receiver. The performance difference can be substantial for small to moderate B. 
With 1 feedback bit per degree of freedom, the MMSE receiver outperforms a matched filter by 
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as much as 30%. However, an MMSE filter is more complex than a matched filter. Therefore, 
there is a performance tradeoff between feedback and receiver complexity. 

We also simulated a multipath fading channel in which each user's signal transverses 2 paths 
with different gains {E\hk^i\'^ = 0.9 and i?|/ifc2p = 0.1, V/c). Furthermore, K interfering users 
are divided into 2 groups. Ki users transmit signal with = Pi while K2 users with = P2- 
This scenario may follow from a system with differentiated quality of service. We obtain the 
large system SINR from Theorem [3] with the asymptotic distribution of Ai 

fA, (a) = ^6ia -P,) + ^5{a- P^) (60) 

where normalized loads Ki = Ki/N and K2 = K2/N. Both the large system and corresponding 
simulated results with Ki = K2 = 0.25 and different sets of Pi and P2 are shown in Fig. IH The 
large system performance closely approximates the performance of the system with = 32. As 
grows, the performance of a finite-size system will converge to that of the large system. In 
this example, reducing the transmit power of one group of users by 20 dB (P2 from 10 to 0.1) 
decreases the required feedback to achieve 0.5 dB away from the single-user performance by 
B = 0.4. 

V. Conclusions 

We have shown expressions for a large system SINR for RVQ with both matched filter and 
linear MMSE receiver. The SINR is a function of normalized load (number of users per degree of 
freedom) and normalized feedback bit (number of feedback bit per degree of freedom). Both ideal 
nonfading channel and multipath fading channel were considered. The SINR of the quantized 
signature for both receivers increases exponentially with B. For a small load, RVQ achieves 
close to the single-user performance with only fraction of feedback bit per quantized signature 
coefficient. We compared performance of the MMSE receiver with that of matched filter derived 
in [9] and showed that the performance gap is large for small B. The simpler matched filter 
requires more feedback to achieve a target SINR than the MMSE receiver does. 

In this work, we assume that the receiver can estimate channel and interference covariances 
perfectly. In practice, a very accurate channel estimation is achieved by a large amount of 
training. How the performance of RVQ is affected by imperfect channel estimate at the receiver 
(or limited training) was studied by [24], [25]. Here we consider signature quantization for a 
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single user. Future work includes performance analysis of group of users with RVQ-quantized 
signatures. 

Appendix 
Proof of Corollary [2] 

We rewrite (|43] ) as follows 

*(P,/?rvq) = / log(l +p(/3,^q - -^))4,st(x) dx (61) 

log(x + al + piP^qix + al) - l))fs^sl (^) ~ / ^^^^^ + ^n)4isl (^) 

(62) 



where 



/ N \/ ix — a)(b — x) 

(x) = for a<x<b, (63) 

^1^1 27rx 



where a = (1 - and 6 = (1 + v^)^ for fsT > 1. 

To determine p*, we take the first derivative of (l62l) with respect to p given by 

d$(p,/?-)_l 1 1 /• 1 



dp P p(P/9rvq + l) (P/5rvq + l)\i X-y'^^^i 



4^.t(x) dx (64) 



where is the Stieltjes transform of /g^gt(-) and 

We solve for p* (or equivalently y*) by setting (|64l) to zero and obtain 

= (p*/3,^q + (66) 

Substituting the Stieltjes transform of g{-) and using the change of variable from (|65l) in (l66l) 
give 



-1 + /sT - ± v/(y*)2 - 2(is: + i)y* + (fc - 1)2 _ 



rvq 



(67) 



2y* l_/5-(^* + a2) 

Simplifying (l67l) gives 

. ^ (l-/3-(ir-l + aD)(l-/3-a^) 

^ /5r^q(l-/3r^q(^+^^)) ' ^ ^ 



With change of variable (I65I) . we obtain 



^ ^ 1 L (69) 

/^rvq(-'- /^rvq'^n) (/^rvq) /^rvq 
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To show that p* achieves the maximum, we compute the second derivative of <^{p,(3^^ (|621 ) 
with respect to p 

d2$(p,/5-) (a; + - 1) 



dp2 



< 0. (71) 



For large enough > p^^*, p* ^ can exceed I / {P^^^ - P^q) ■ To determine P^^*, we set 

K 111 

(72) 



Simplifying (l72l) gives the following quadratic equation 



[K + al- Pa'jP^^' + [(1 -K- a'JP - 1]P^^ + P^. = 0. (73) 



Solving (|73]) gives the only solution 

{K - VWy + 2al{K - Vk) + < 

Thus, 



/?rvq* = 7- : " :• (74) 




P* = / MrvqVJ^-Mrvq^n; VMrvq;" Mrvq ''^'^ ^^rvq (75) 

Substituting p = p* in (1621) and rearranging give 

$(p*,/3- ) = log(p*/3,^q + 1) + /'log(x + - f )4^^t(3:) dx 



rvq 

2\ 



(76) 



log(x + a„)/g^gt(x) dx. 

First, we consider the case where P < P^^ < P^^*. Substituting p* into the first term in (|76l ) 
gives 

log(p*/?.vq + 1) = log(T3^ - (77) 

To evaluate the two integrals in (1761) . we apply the following lemma. 
Lemma 1 ( [26]): For K > 1, 

^ log(x + a)4^st(a;) dx = log(w;(a)) - {K - 1) log(l - ^v{a)) - v{a) (78) 
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where 



w{a) = ^{l + K + a+^^/{l + K + ay- AK), (79) 



v{a) = -{l + K + a-^{l + K + a)^-AK). (80) 
Using Lemma [U and (|77]), we can evaluate (|76l) for (3 < [3^ < (3^ * 



where 



P = T>7^ , ^ T - 7^ ^cr„. (82) 



1 - /^rvq^n 1 
- 1 + a2 /?, 

Next we evaluate $(p*, /^^^^^q) for 13"^^ > P^^*. Substituting the value of p* from dVS]) gives 

log(p*/3,^q + 1) = log(/3:S.x) - log(/5:S.x - /5rvq) (83) 

and 



7 = - :5 = -(1 - yKy. (84) 

r A^rvq A^max 

Substituting (l84l)into the second term in (1761) and applying Lemma [T] gives 
/'log(x + al- -^_^)4^^t(x) dx = ^KlogiK) - {K - 1) log(v^- 1) - v^. (85) 



Thus, for [3^^ > (3^ * 



rvq 

rvq ^ A^rvq 



$(p*, ) = log(/?^,J - log(/3~ . - P^^) + log(ir) - (ir - 1) log(v^ - 1) - 

- log(«;(a2)) + (ir - 1) log(l - ^v{al)) + t;(a^). 



(86) 



Also <l>(p*,/?r^q) = i?log(2). We can explicitly solve for [3^^ as follows 

/?rvq = /3max(l - 2-^[exp{liriog(ir) - (i? - 1) M ^"^"^ ) 

-log(^(cr^))+^;(o^D-V^}]). 
To solve 5* which is corresponding to (3"^^ dUl), we substitute (3"^^ in (l87l) . 



(87) 
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SNR = 5dB 




Fig. 1. Shown is a large system SINR for MMSE receiver versus normalized feedback bit B with different normalized loads 
K = 0.25, 0.5, 1, 1.25 and SNR = 5 dB. 
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K/N = 0.5;SNR = 8dB 




2 





Interfering signatures contain i.i.d Gaussian elements. 

' - ' - ' Interfering signatures are ortliogonal. 




1 1 1 1 1 



0.5 1 1.5 2 2.5 3 

B/N 



Fig. 2. Shown are large system SINR's for different distributions of interfering signatures with K = 0.5 and SNR = 8 dB. 
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Fig. 3. The large system SINR for MMSE receiver is compared with the approximation derived in [8] and the large system 
SINR for matched filter [9]. Also shown is the simulation result for N ^ 12, K = 0.75 and SNR = 10 dB. 
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Fig. 4. A large system SINR for MMSE receiver and multipath fading with two groups of users is shown with simulation 
results. SNR = 5dB, number of paths L = 2 for all users, and K = 0.5. 
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